Dusty Disks around White Dwarfs I: Origin of Debris Disks 



Ruobing Dong^ ^*, Yan Wang^ ^ D. N.C. Lin' \ and X.-W. Liu^-^ 

^Kavli Institute of Astronomy and Astrophysics, Peking University, Beijing, China 

^Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA; 

rdong @ astro.princeton. edu 

^Department of Astronomy and Astrophysics , Penn State University, State College, PA 16802, 

USA; yuwl23@psu.edu 

'^Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064, 

USA; lin@ucolick.org 

^ Department of Astronomy, Peking University, Beijing China; liuxw@bac.pku.edu.cn 

ABSTRACT 

A significant fraction of the mature FGK stars have cool dusty disks at least an 
orders of magnitudes brighter than the solar system's outer zodiacal light. Since such 
dusts must be continually replenished, they are generally assumed to be the coUisional 
fragments of residual planetesimals analogous to the Kuiper Belt objects. At least 10% 
of solar type stars also bear gas giant planets. The fraction of stars with known gas 
giants or detectable debris disks (or both) appears to increase with the stellar mass. 
Here, we examine the dynamical evolution of systems of long-period gas giant planets 
and residual planetesimals as their host stars evolve off" the main sequence, lose mass, 
and form planetary nebula around remnant white dwarf cores. The orbits of distant gas 
giant planets and super-km-size planetesimals expand adiabatically. During the most 
intense AGB mass loss phase, sub-meter- size particles migrate toward their host stars 
due to the strong hydrodynamical drag by the intense stellar wind. Along their mi- 
gration paths, gas giant planets capture and sweep up sub-km-size planetesimals onto 
their mean-motion resonances. These planetesimals also acquire modest eccentricities 
which are determined by the mass of the perturbing planets, the rate and speed of stel- 
lar mass loss. The swept-up planetesimals undergo disruptive collisions which lead 
to the production of grains with an extended size range. The radiation drag on these 
particles is ineff"ective against the planets' resonant barrier and they form 30-to-150- 
AU-sizes rings which can eff"ective reprocess the stellar irradiation in the form of FIR 
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continuum. We identify the recently discovered dust ring around the white dwarf WD 
2226-210 at the center of the Helix nebula as a prototype of such disks and suggest 
such rings may be common. 

Subject headings: stars: AGB and post-AGB — stars: evolution — (ISM:) planetary 
nebulae: general — stars: individual (WD2226-210) — (stars:) planetary systems — 
methods: numerical 



1. Introduction 

The discovery of over 300 planets shows that at least 10 % of nearby solar type stars have 
Jupiter-mass planets around them (Gumming et al. 2008). According to widely adopted sequential 
accretion hypothesis, these planets formed through condensation of heavy elements, emergence 
of planetesimals which coagulated into protoplanetary embryos and cores prior to efficient gas 
accretion. Through a series of population-synthesis simulations based on this scenario (Ida & Lin 
2004a,b, 2005, 2008a,b), we were able to reproduce the observed mass-period (Mp-P) distribution 
of the known gas giant planets (Schlaufman et al. 2008) around a modest fraction (t]j) of FGK main 
sequence dwarf stars. 

Mass accretion rate M from protostellar disks onto stars more massive than the Sun is ob- 
served to increase with square of the stellar mass i.e (Natta et al. 2006, Garcia-Lopez et al. 
2006, cf Clarke & Pringle 2006). If the initial mass of these disks also increases with M*, the 
greater availability of heavy elements would provide more rapid growth and larger of protoplane- 
tary embryos around intermediate mass stars (Ida & Lin 2004b, 2005, Kennedy et al. 2007). Fur- 
thermore, the extended local surface density and pressure maxima at the ionization radius (where 
the mid-plane temperature - 1000 K), snow line (Kretke et al. 2009), and outer edge of the 
dead zone provide more efficient barrier against hydrodynamic drag for the sub-km grains and 
type-I migration for the critical-mass (~ lOM^) cores. The enhance retention efficiency of these 
building blocks enhance emergence probability of gas giants (Laughlin et al. 2004, Ida & Lin 
2004b). 

On the observational side, high-precision radial-velocity surveys of intermediate-mass main 
sequence stars have been handicapped by their hot and active atmosphere as well as their fast spin 
(Griffin et al. 2000). However, these problems are significantly reduced when these stars evolve 
off their main sequence. Several planets have been found around G sub-giants and giants (Sato 
et al. 2005) which have more expanded envelopes and relatively cool atmospheres. Preliminary 
surveys indicate that the fraction of these relatively massive G giant stars with gas giant planets 
('/z ~ 1 /3) may be more than double that around the solar-type G dwarf main sequence stars (Sato 
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et al. 2007, Johnson et al. 2007, DoUinger et al. 2007, Hatzes 2008, Johnaon 1008). 

Another area where observations provide valuable clues on the planet formation process is the 
search for residual building-block planetesimals around their host stars. Since 1995 a population 
of Kuiper Belt objects (KBO's) has been found in the solar system and beyond the orbit of Neptune 
(Jewitt & Luu 1995, Brown et al. 2005). These objects are thought to be the parent bodies of a 
diffuse dusty ring which emits cold zodiacal light with far infrared (FIR) radiation. 

A series of recent systematic 70 ^ira (FIR) observational surveys with the Spitzer Infrared 
Space Telescope indicate that 1) none of the nearby M stars show detectable excess (Gautier et al. 
2007), 2) 15 % of 274 FGK stars show excess (Bryden et al. 2006, Beichman et al. 2006, Trilling 
et al. 2008), 3) 30 % of 160 A stars also show excess (Rieke et al. 2005; Su & Rieke 2007). The 
threshold level for a marginal detection by the Spitzer Telescope is at least an order of magnitude 
more intense that of the solar system zodiacal light. 

Since both the fraction of stars with gas giants (rjj) and that with detectable debris disks (t/j) 
increases with the stellar mass (M*), we can extrapolate that a major fraction of intermediate- 
mass stars may have both gas giants and persistent, rich debris disks. This inferred association of 
gas giant planets with debris disks is highlighted by the recent discovery (with adaptive imaging 
technique) of three gas giants and a debris disk outside them around an A5V star HR 8799 (Marois 
et al. 2008). The mass and periods of these three planets are (10, 10, 7) Mj and (24, 38, 68) AU. 
Around this = 1.5Mq star, the ice line is located at ~10 AU. In §2, we provide a more detailed 
discussion on the association of debris disks and planet formation in the Solar System and around 
nearby stars. 

With the evidence of existence of many detectable debris disks and gas giant planets around 
nearby, mature, intermediate-mass stars, we now consider their dynamical evolution during the 
post main sequence evolution of their host stars. The goal of our study here is to identify whether 
the inconspicuous planet building blocks may reveal their presence. The analysis to be presented 
below is particularly relevant for the AGB mass loss phase where the mass loss rates range up to 
~ 10 "^Mq yr"^ and where most of the stars may reach this mass loss rate near the end of the AGB 
phase (Willson 2000). The typical outflow speed is ~ 10 km s ' and may vary on the time scale 
of 10^""* yr (Weidemann 1987, Blocker 1997). After losing mass on the AGB, these stars rapidly 
evolve across the HR diagram to form planetary nebula around emerging white dwarfs. During the 
post AGB phase, although high-speed mass loss continue to occur, the mass flux is much reduced 
(cfKwokl982). 

In this paper, we mostly consider stars with a main-sequence masses in the range of 2Mq < 
M* < 4M© (A-F stars). We choose these stars not only because many, perhaps most of them, have 
gas giants and debris disks (similar to those around the less massive Sun and HR 8799), but they are 
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also thought to be the progenitors of some well known planetary nebulae such as the Helix nebula 
where the mass of the remnant white dwarf WD 2226-210 is ~ O.SSMq (Su et al. 2007). This rapid 
and large-magnitude change in the gravitational potential can lead to great orbital expansions for 
their planets. This evolutionary tendency may lead to the possibility of resolvable HST imaging of 
long-period planets which are otherwise undetectable. 

Evidence for the existence of residual planetesimals comes from the recent discovery of a 
dusty ring around a young white dwarf WD 2226-210 at the center of the Helix nebula (Su et al. 
2007). This ring extends between about 35 and 150 AU from a central white dwarf WD 2226- 
210 (well interior to the complex, helix structure), and have a total mass of about 0.13 M®. We 
propose that this ring is the byproduct of planets and planetesimals' orbital evolution during the 
epoch when the central star rapidly lost most of its mass. Based on the extrapolation of the system 
around HR 8799, our basic conjectures are: 

1) gas giant planets and debris disks are common around A main sequence stars; 

2) at the epoch of their main sequence turnoff, these relatively massive stars have debris disks with 
more than lOM^ in total mass and most of it is contained in km-size residual planetesimals 

3) gas giant planets' orbits expanded adiabatically by nearly an order of magnitude while there 
host stars evolve through AGB mass loss phase with rapid stellar mass loss; 

4) the KBO-equivalent planetesimals are captured onto the main motion resonance of gas giants; 
and 

5) through collisions, these parent bodies generate dust which reprocess the radiation of their 
central star. 

These basic assumptions are based on their solar system and HR 8799 analogues. A detailed 
and systematic analysis of these processes is useful for understanding the long-term evolution of 
this class of planetary systems. In §3, we first construct a working model for a solar-system like 
configuration. The orbital evolution of the planets and planetesimals is computed and analyzed in 
§4. Here we take into account of the central star's mass loss, the hydrodynamic drag by the expand- 
ing envelope and the planet's dynamical perturbation on the residual planetesimals. In §5, we show 
that collisions between the resonant planetesimals lead to their disruptions and fragmentation. We 
derive the condition for the retention of dust grains into rings outside the orbits of some gas giant 
planets. We model the structure of the disk around WD 2226-210 in §6. Finally, we summarize 
our results and discuss their implication, especially in the context of possible detectable gas giant 
planets around young white dwarfs in §7. 
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2. Association of debris disks and gas giants 

In this section, we briefly recapitulate existing observational information on the inventory 
of residual planetesimals in the outer solar system and around mature, intermediate-mass main- 
sequence stars. These data will be used in the subsequent sections for the development of our 
working model. 



2.1. Inventory of residual KBO's in the outer solar system. 

KBO's are marginally preserved at their remote location because planetary perturbation (by 
gas and ice giants) is relatively weak (Duncan & Quinn 1993). However, their present-day spa- 
cial distribution also poses a challenge to the formation of the largest KBO's which have masses 
comparable to Pluto. If their building blocks were initially smoothly distributed with respect to 
their distances from the Sun, in situ formation of such large KBO's would require more than a few 
of planetesimals in the region 30-50 AUs (Stem 1996, Youdin & Chiang 2004). This required 
total mass is ~ 2 orders of magnitude larger than the total reservoir of all known KBO's. Dynam- 
ical perturbation by Uranus and Neptune can induce a substantial depletion from a massive initial 
inventory to the present-day KBO reservoir (Duncan & Quinn 1993). This stringent requirement 
would be relaxed under the assumption that the most massive KBO's were formed closer to the 
Sun and then scattered to their present-day location by gas or ice giant planets either during the 
advanced phase of their formation (Zhou & Lin 2007) or through long term dynamical evolution 
of the system (Duncan et al. 1987). 

It is possible that the size distribution of the KBO's extend well below the present detection 
limit. Elsewhere in the solar system, the crater density on the unprocessed moons of the outer 
planets (Hartmann 1999) indicate that they have been bombarded by meteorites with a size distri- 
bution dN/ds oc s~^-^ as observed among the asteroids (Gradie et al. 1989). This distribution is 
probably the byproducts of a coUisional equilibrium (Donanyi 1969, Williams & Wetherill 1994, 
Farinella & Davis 1996, Tanaka & Ida 1996, Kenyon & Bromley 2002, Weidenschilling 2004). 
On the small-size limit, particles in Saturn's rings (Cuzzi et al. 1980) and grains in the interstellar 
medium (Mathis et al. 1911, Laor & Draine 1993) follow a similar pattern. With such a size 
distribution, most of the surface area is covered by the small particles while a large fraction of the 
total mass is contained in the large particles. 

As a population, the KBO's have a directly-observed size distribution dN/ds oc s"^ for s > 
100km (Trujillo et al. 2001) whereas, below .y ~ 70 km, dN/ds is considerably shallower than 
s~^-^, (i.e. there is a deficit of small particles) down to the detection limit of a few 10 km (Bernstein 
et al. 2004). This break in the slope of the size distribution implies a large fraction of the KBO's 
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total mass is contained in particles with ~ 10 - 100 km. One possible reason for this complex 
power-law distribution is the coUisional fragmentation of KBO's as piles of loose gravels (Pen & 
Sari 2005). 

At their present distance, it is difficult to directly image planetesimals smaller than a few 
10km. However, micron to sub-mm particles in the asteroid belt and Kuiper belt scatter and repro- 
cess the solar radiation to produce warm and cold zodiacal light (Backman et al. 1995). Due to the 
Poynting-Robertson (PR) drag induced on them by the stellar photons, the grains undergo orbital 
decay on time scale much shorter than the age of the Sun. Thus, these grains must be continually 
replenished, most likely through the break up of their parent bodies in the Kuiper Belt (Kenyon 
& Luu 1999). Energetic collisions can indeed lead to the disintegration of parent bodies and a 
catastrophic production of very small fragments (Fujiwara et al. 1911, Petit & Farinella 1993, 
Benz & Asphaug 1999). The KBO's can provide sufficiently frequent collisions to replenish the 
supply of small grains needed to account for the observed level of the zodiacal light provided 1) 
a non-negligible fraction of their total mass is contained in KBO's with < 1 - 100 km, and 2) a 
modest fraction of these collisions actually lead to catastrophic fragmentation. 

Asymmetric reradiation of stellar light by spinning bodies can lead to a torque on their orbits. 
This effect is commonly known as the Yorkovski (YORP) effect (Bottke et al. 2000, de Pater & 
Lissauer 2001). The characteristic YORP-induced orbital-evolution timescale for planetesimals 
with a density Pp, semi major axis a,, and Keplerian speed vk = {GM^/ajy^^ is 

Ty ~ ncppSa^VK/fL, (1) 

where M, and L* are the mass and luminosity of the host star. The magnitude of the efficiency 
factor / is determined by the fractional amplitude of seasonal and diurnal temperature variations. 
The orbits of planetesimals with prograde spins expand as they experience a positive force (and /) 
whereas planetesimals with retrograde spins undergo orbital decay with a negative /. The YORP 
effect enhances the PR drag for small particles. It can also lead to significant diffusion in the 
spacial distribution of km-size planetesimals during the main sequence life span of intermediate- 
mass stars. At 10 AU from stars with L, ten times larger than that of the Sun, meter-size particles 
may undergo orbital decay within a Myr (see O below). 



2.2. Debris disks around nearby intermediate-mass stars 

In the quest to determine the properties and understand the origin of planets around other 
stars, it is often useful to identify their solar system analogs. Although it is impossible to directly 
detect residual planetesimals similar to the KBO's, the counterpart to the solar system zodiacal 
light have been found with Spitzer Space Telescope in the form or MIR and FIR excesses around 
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a significant fraction of nearby mature solar type stars (Rieke et al. 2004, Beichman et al. 2005, 
Moro-Martin et al. 2007, Hillenbrandt et al. 2008). Such excesses are assumed to be due to the 
reprocessed stellar radiation by debris dusty disks or rings at ~ 2 - 20 AU from their host stars 
(Bryden et al. 2006, 2007, Trilling et al. 2008). The relevant-size grains are either blown away 
by radiation pressure or spiral toward their host stars due to radiation drag (both PR and YORP 
effects) on time scale much shorter than the age of the system (Artymowicz & Clampin 1997). 
They must be continually replenished by the collisions between relic planetesimals, analogous to 
the those in the Kuiper Belt and Oort Cloud (Kenyon & Bromly 2002, 2004, Wyatt et al. 2007, 
Moro-Martin et al. 2007, 2008). 

The Spitzer Space Telescope's detection threshold for these system is at levels at least an 
order of magnitude brighter than the outer solar Zodiacal light. It is entirely possible that all stars 
may have debris disks similar to that around the Sun. Some known debris disks have FIR excesses 
are 10^"^ times larger than that associated with the zodiacal light in the solar system (Bryden et 
al. 2006). On average, the rate of coUisional dust production is approximately proportional to the 
square of the surface density of their parent bodies. Presumably, stars with detected debris disks 
have richer reservoirs of relic planetesimals than the sun, albeit occasional major impacts may lead 
to brief episodic flare ups in the IR excess around individual systems (Kenyon & Bromly 2004, 
Grigorieva a/. 2007). 

In the attempt to differentiate transient events and average properties, it may be fruitful to 
observationally determine the fraction of stars with detectable debris disks in a given population- 
analysis sample. In §1, we have already cited observational indications that rjd is an increasing 
function of M*. This correlation of debris disks with the stellar mass can be partially attributed to 
observational selection effect because 1) A main sequence stars are generally younger and 2) their 
disks with similar fractional luminosity, relative to theirs, are brighter, warmer, and easier to detect 
than those around FGK dwarfs (Bryden, private communication). Nevertheless, this correlation is 
consistent with the expectation inferred from 1) the fraction of stars with planets, and 2) the disk 
accretion rate and mass around 1-3 Myr old stars are increasing functions of the stellar mass. 

Further evidences for the coUisional origin of debris-disk dust particles can be inferred from 
the general decline in the intensity of the debris disks with the age of individual host stars (Rieke et 
al. 2004). As a population, the fraction of A stars with 70 yum excess decreases from 48 % around 
those younger than 30 Myr to 12 % for those older than 400 Myr. These signatures are consistent 
with the expectation that both dynamical and coUisional processes lead to the depletion of these 
parent bodies. Nevertheless, a fraction of relatively massive A stars returns detectable debris disks 
throughout their main sequence lifespan. 

The fraction planet-bearing FGK main sequence stars with FIR signature appears to be com- 
parable to those without any known Jupiter-mass planets. After correcting some unfavorable ob- 
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servational selection effects, the intensity of the FIR excess for gas giant planet-bearing stars is ~ 3 
times that for stars without any planets, albeit the statistical significance of this correlation is weak 
(Bryden et al. 2007). In contrast to the strong correlation between the metallicity of FGK dwarfs 
and their planet-bearing probability (Fischer & Valenti 2005, Santos et al. 2005), there appears 
to be no correlation between it and detectable debris disks (Greaves et al. 2007). This dichotomy 
is probably due to the retention of debris-disk parent bodies depends much less sensitively on the 
surface density of planet-building blocks in their nascent disks (Wyatt et al. 2007) than the emer- 
gence of gas giant planets (Ida & Lin 2008a, b). However, the relatively high fraction of stars with 
planets around more massive (with mass M* > 1 - 2M©) stars is independent of their metallicity 
(Pasquini et al. 2007). We have already indicated above that planet formation around the relatively 
massive stars may be more prolific not only because the protostellar disks around them have greater 
inventories of planet building blocks but also that the planetesimals formed more readily and are 
better retained in these disks. 

Since both rjj and rjd increases with M, and attain modest values (0.2 - 0.3) for intermediate- 
mass stars, a significant fraction of them (such as HR 8799 and Fomalhaut Kalas et al. 2008) may 
contain both. Our objective in this paper is to consider the dynamics of the planets and residual 
planetesimals during the post main sequence evolution of some intermediate-mass host stars. 



3. A working model 

Dynamical response of planets as a consequence of mass loss during the planetary nebula 
stage has been considered by several authors in the past mostly in the context of metal-rich DZ 
white dwarfs (Stem et al. 1990, Sackmann et al. 1993, Parriott & Alcock 1998, Siess & Livio 
1999a, b, Debes & Sigurdsson 2002, Jura 2006). 

For the origin of dusty disks in planetary nebulae, we consider similar effects, with a greater 
detailed treatment of both orbital dynamics and gas drag effects. We idealize the dynamics of the 
system into a series of interaction between three populations of objects: planetesimals with mass 
ruj, planets with mass Mj, and a central star with mass M*. 

There are two mean factors which influence the orbital change of a planet around the central 
star during the epoch of its mass loss: the change of the star's mass, which reduces its gravity and 
causes the planet's orbit to expand; and the gas drag by the outflow, which induces orbital decay. 
But the magnitude of this drag on a gas giant planet is negligible so that the equation of motion for 
the gas giant is reduced to the standard Kepler's problem. 

But the drag effect has a significant effect on the orbit of the planetesimal. In addition, it is 
also perturbed by the gravity of the planet. In a non rotating frame centered on the star, the equation 
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of motion of a planetesimal at a position r,- can be expressed as 

Svt _ GM.r, ^ GM/r, -r,) ^ GM,r, f, 
df- ■' |r,--rjp ■' r'. nii 

where is the position of the planet and M, is a function of time. On the right hand side of the 
above equation, gravity of the star and planet are expressed in the first two terms. The third term 
corresponds to the indirect force due to the rotation of the frame (Murray & Dermott 1999) and the 
forth term represents the gas drag force on the particles. 

In the Stoke's region, the magnitude of the drag force on a planetesimal by the wind of its 
host star can be expressed as 

= Wwindis^ + sl){ri - y„ind)\ri - y„ind\ (3) 

where v^nd is the wind velocity, s = (Smi/Anppy^ and Sg = Gmj/v^^ are the physical and gravita- 
tional radius of a planetesimal with a mass m and density pp (cf Supulver & Lin 2000). In the radial 
direction, this drag contribution is much smaller than the gravity from all three bodies. However, 
in the azimuthal direction, the drag slows down the planetesimals' initial Keplerian orbits (with a 
velocity Vk = (GM^/ry^^ at a distance of r) to migrate inward on a characteristic time scale 



^drag 



Vk 



2vk 



2fi ~ 3M, Vk ' 



In the above expression, we have neglected the gravitational focussing eff"ect (which is small 
for km-size planetesimals) and assumed a steady spherically symmetric wind such that the back- 
ground gas is determined by the stellar mass loss rate 

Mh. = 4-nr^p^indVmnd (5) 

and the wind speed Vmn^i is generally assumed to be lOkm/s which is larger than vk for r greater 
than a few AU. 

We can include other effects to the equations of motion. For example, the shape of some 
planetary nebulae significantly departs from spherical symmetry (O'Dell et al. 2004). In a follow 
up paper, we add a J2 component to the stellar gravity. It is also possible to include the gravity 
of additional stellar or planetary companions. In the current investigation, we neglect the PR and 
YORP effects during the AGB mass loss phase. Based on eq©, we estimate that these effects can 
lead to significant orbital evolution during the brief AGB mass-loss phase (for ~ lO'^ yr) only for 
sub-meter-size particles. Comparing eqs ([T]) and (Hj), we find that the hydrodynamic drag has a 
stronger impact on the orbital evolution of all particles at all locations if 

M>My = %fUv^,ndlcv\ - lO-^Meyr-i (6) 
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during the AGB mass loss phase. 

However, during the main sequence life span t» of A stars, YORP can induce residual plan- 
etesimals with sizes smaller than 

Sy = TjLJ{ncppa]vK) ~ {aillQAUf'^km (7) 

to migrate over a significant fraction of their a,. The inward migration of planetesimals may be 
halted after they are captured by the mean motion resonances of any embedded gas giants whereas 
the outward migration is limited to be less than 10^ AU (see further discussions on the initial 
planetesimal reservoir below). 

In order to fully explore possible range of initial conditions, we solve the equation of motion 
with a Hermit scheme which is kindly provided by Dr Sverre Aarseth (2003). 

4. Orbital evolution of planets and planetesimals due stellar mass loss 

We recapitulate here the response of planets and planetesimals due to the changing gravita- 
tional potential during the rapid loss phase of their host AGB or post- AGB stars. 



4.1. Adiabatic orbital evolution of planets and large planetesimals 

In our model, there are two main factors which aff'ects the planet's orbital evolution: the 
reduction of the central mass and the gas resistance. The resistance term increases with while 
the gravitational term increases with m, ~ . For relatively large planetesimals, gas resistance 
is unimportant and its orbital change is determined only by the central mass. Provided the time 
scale of the host star's mass loss Tm = |M,/M*| is long compared with their Keplerian orbital time 
scale tk = {a^GM^y^, the planetesimals' angular momentum h = vR does not change during this 
episode of mass loss, i.e. it is an adiabatic invariant. The planetesimal also retain their radial action 
so that a centrifugal balance is maintained in which = GM*/fl,. Consequently, a, = h^/GM and 
the orbit expands as the host star loses mass (Duncan & Lissauer 1998). 

In contrast, distant planets and planetesimals with tk > Tm react to the impulsive change in 
the gravitational potential such that they would escape if their host stars lose more than half of their 
mass. During the AGB mass loss phase, the characteristic time scale for mass loss, tm ~ 10'^~^ yr. 
All planets and planetesimals with period less than tm (or semi major axis within several 100 AU) 
undergo orbital expansion rather than escape. 
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The characteristic orbital expansion time scale is 



' exp 



(8) 



From equations (HI) and dH), we find that the magnitude of Texp - Tdrag for particles with sizes 



3M, 



Vk 



^ 0.1 — 



1.5 



Ig cm 
Pp 



lOAU 



^wind 



km. 



(9) 



The above equation is applicable in the p„,,>,rf oc 1 /r^ (or v„ind = constant) region of the stellar wind, 
and with the planet orbital radius greater than a few stellar radii. During the AGB mass loss phase, 
orbits of planetesimals with s > Sc expand during the stellar mass loss. Orbits of planetesimals 
with s « Sc undergo orbital decay on a time scale 
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The relation of Sc with orbital radius is showed in Figure [U which indicates a km-size planetes- 
imal would be larger than 5c at r = 15 AU from a mass-losing star so that its orbit would expands 
whereas it would be smaller than at r = 10 AU so that it would undergo orbital decay during the 
mass loss. This dependence arises because the gas density in a steady wind decreases with a,-. All 
planetesimals with sizes larger than Sc at their initial location are expected to expand, albeit with 
size-dependent speeds. 

We are primarily interested in Kuiper-Belt equivalent regions because 1) that is the region 
where FIR from the debris disks implies the potential presence of a large population of residual 
planetesimals, 2) it is outside the region where gas giants are expected to form, and 3) it is also 
outside the quasi hydrostatic atmosphere of AGB's so that planetesimals may survive the evolution 
of their host stars prior to the AGB mass loss phase. For this region, Sc ~ I km. In §2, we have 
already indicated that a significant fraction of the total mass in the Kuiper Belt may be contained 
in KBO's in this size range. In the above section, we also indicated that due to the YORP effect 
during the mean sequence evolution of their central stars, km-size KBO's (with s < sy) have a 
tendency congregate near the outer mean motion resonances of long-period gas giants. 

Around any given host star, the mass loss rate during the AGB mass loss phase may have large 
amplitude fluctuations. Since both the orbit expansion in response to the decreasing stellar mass 
and the gas drag in the stellar wind are proportional to the stellar mass loss rate, Sc is independent 
of M [note that Sc is defined by Texp - "^drag and both Texp and T^rag are inversely proportional to M]. 
Nevertheless, the magnitude of Sc decreases with the stellar mass. Also, tk is large at relatively 
large a, such that the condition for adiabatically may be violated during episodes of rapid mass loss 
when tm is reduced below tk- Nevertheless, planets and planetesimals can still be retained at large 
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distances provided the fractional loss of mass is limited during these brief events. Substantial mass 
reduction of the host stars can also cause dynamical instability of multiple-planet systems (Debes 
& Sigurdsson 2002). We will further explore this possibility in a subsequent paper. 



4.2. Entrainment of small planetesimals by the gas outflow 

We demonstrate the orbital evolution of marginal and small planetesimals with numerical 
solutions of the equation of motion. In general, the orbits of intermediate (i' ~ 1 km) size, more 
distant (r ~ 10 AU) planetesimals may enlarge at rates slower than the adiabatic expansion rates 
for the more massive and distant planetesimals. 

If a planetesimal is comparable to or smaller than a critical size, orbital decay would be 
introduced due to gas drag. The gas drag would reduce the planetesimal 's momentum notably, and 
the angular momentum would not remain constant. From equations ^ and (O, we find that drag 
by the outwardly flowing gas exceeds the stellar gravity on particles smaller than 

2 

^ ^PwindV„ind M^Pk Vwind 
■^blowout ~ ^ ~ ^ V ^ / 

where Pk is the orbital period. These very small particles {s < i'biowout ~ 1 cm) are entrained by 
the stellar wind. 

In order to quantitatively illustrate the divergent evolutionary paths of these particles, we 
consider a series of 4 planetesimals with radius of 10m, 100m, 1km and 10km all starting at 15Au 
from the star, whose orbits are initially without eccentricity. The mass of the star changes linearly 
from an initial value M, = 3 Mq to a final value Mf = 0.5 Mq in 50,000 years. The orbital evolution 
of these planetesimals is shown in Figure [2l The 10 km object's has a final orbital radius 6 time 
of the its initial value, which is exactly the reduction factor of the stellar mass. This correlation 
implies that the planetesimal's orbital evolution is adiabatic. The 5=1 km object also moves 
outward, but its expansion factor is considerably smaller than that of the 10 km-size planetesimal. 
The 100m object undergoes a small amount of orbital decay and its orbital evolution is stalled at 10 
AU. The 10 m planetesimal quickly spiral into the central region and it is expected to be evaporated 
at about 10,000 years. 

The processes of orbital decay (for planetesimals with 5'biowout < s < Sc) and entrained outflow 
(for planetesimals with s < ^blowout, due to the hydrodynamic drag by the outflowing gas, are 
equivalent to the Poynting-Robertson (PR) drag and radiative blow out by stellar radiation. During 
the AGB mass loss phases, hydrodynamic efl'ects are much more intense than radiative efl'ects. 



-13- 



4.3. Capture of small planetesimals by the planets' mean-motion resonances 

In the above example, we neglected the presence of any major planets. Such planets also 
undergo adiabatic expansion during the planetary- nebula phase when their host stars lose a sub- 
stantial fraction of their masses. For these large entities, the orbital expansion is adiabatic and take 
place on a time scale r^xp - r^. In addition, they exert dynamical perturbations on the nearby plan- 
etesimals, especially those on their low-order mean-motion resonances where their orbital periods 
are commensurate with each other in terms of small integers. 

During the stellar mass loss, planets' orbits expand from their initial location to their final 
location r^nai = (Mi/Mf)rp. Beyond these planets' orbits, small or marginal-size planetesimals, 
initially located at a, (we will stick to r for the radius of planet, since it always keeps a circular 
orbit, and a for the semi-major axis for planetesimals.), either undergo orbital decay or expand at 
slower paces. Since their differential motion is determined by the gas drag effect, a planet and a 
(more distant) planetesimal cross each other's orbit on a time scale 

(^i - fp) (ar-rp)s 

Ta - Tdrag(^) = Tm- (12) 

Clj CI J S 

Provided their initial location a, < rfinai, particles with sizes between 

Scat ^ 

■^catchup ~ ^ N vJ^'^/ 

(«/ - rp) 

and ^blowout generally undergo orbit crossing with the planet because the expansion of their orbits 
is slower than that of the planet. 

Planetesimals along a planet's expansion path encounter its mean-motion resonances prior 
to their orbit crossing. Provided the duration across the width of these resonance (A) is long 
compared with the libration period near their centers (rub) planetesimals would be captured onto 
the planets' mean-motion resonances (Murray & Dermott 1999). For example, a planet's 2:1 
resonance has a finite effective width A/rp ^ 2expl/2(A/p/A/*)'^^, where Pk = 2nTK is the orbital 
period of the planet, and a libration period at its center Tiib - (M,/Mp)^^^/'/f/2exp(l/2). For this 
2: 1 resonance, the critical condition for capturing a planetesimal (with a size s) is 

Mp > - . , M,. (14) 
s 4TMexpl 

Planetesimals with ^catchup < * < ^'^scap are captured onto the planet's 2:1 resonance where 

4TMexpl \.6rp Mp 

In the above expression, the term aj/l.6rp is introduced to take into account of the stellar mass 
loss which is the cause of planet's orbital expansion. Equation (fT4l) indicates that prior to reaching 
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~ 50 - 80 AU, a Jupiter-mass planet would be able to capture exterior km- size planetesimals 
(s - i'c(~ 1 km)) onto its 2:1 resonance. An Neptune-mass ice giant would also be able to capture 
these marginal-size planetesimals while its semi-major axis is < lOAU. 

Planetesimals with initial semi-major axis a, < 1.6rp are already interior to the planet's outer 
2:1 resonance and they cannot be captured onto it. Also, relatively small (s < Sc) planetesimals are 
unlikely to be captured onto the planet's 2: 1 mean motion resonances because their orbital decay is 
too rapid. However, this special location is the most distant lowest-order mean-motion resonance 
from the planet. There are many other lowest-order (such as 3:2 and 4:3) mean-motion resonances 
which are closer to the planet's orbit. At these location, planet's resonant capture probability is 
enhanced. 

In the proximity of the planet, these lowest-order mean-motion resonances overlap with each 
other and planetesimals with negligible eccentricities and A ^ (12)^^^i?« can enter into the planet's 
Roche radius (Rr = {Mpl^M^y^^rp) within a synodic period Tsyn - (3rp/2A)PK (Lissauer 1987). 
Using a similar argument as above, the capture criterion becomes ta > Tsyn, i-^- the planetesimals 
must not be able to diffuse through this "feeding zone" before they had an inferior conjunction. 
Thus, a fraction of planetesimals with 



would be able to bypass the planet's expanding orbit. For Jupiter-mass planets at ~ 30AU around 
typical intermediate-mass AGB stars, ^bypass ~ (1 - 10)5^- 

Although YORP effects can modify the radial distribution of km-size planetesimals during 
the host stars' main-sequence phase (see O), we neglect it in the determination of ^catchup, ■^rescap, 
and i'bypass- This approximation is adequate for the AGB mass loss phase provided the host stars' 
M > My (see eq[6l). During the mass-loss process, a planet's orbit expands by a factor of Mi/Mf. 
Planetesimals with sizes in the range 



are shepherded by the planet along its path. This size range is approximately an order of magnitude 
larger and smaller than Sc. A population of slightly larger (^rescap < s < ^catchup) and more distant 
(with initial a, > 1.6rp) planetesimals are captured onto the planet's 2:1 mean-motion resonances. 
Smaller or closer planetesimals may be also captured onto its other (such as 3:2 and 4:3) mean- 
motion resonances. 




(16) 




(17) 
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4.4. Numerical calculation of resonant capture process 

In order to illustrate the resonant capture process, we carried out a series of numerical calcu- 
lations. In Figure [31 we plot the orbital evolution of a 0.1 km-size planetesimal. It has a density 
Pp = 3 g cm"^ and it starts out with a circular orbit at 35 AU around a M, = 4Mq star which linearly 
reduces its mass to Mf = \Mq on a time scale of 100,000 yr through a wind with a speed v„i„d = 10 
km s ^ A planet with a mass Mp = 1 Mj was placed in a circular orbit at 20 AU initially. We 
also present the orbital evolution of the planetesimal when there is no giant planet in the system 
for comparison. 

As mentioned before, the orbital radius of the giant planet only depends on the central stellar 
mass and initial orbit and adiabatically evolve from 20 AU to 80 AU. The planetesimal's semi- 
major axis first decreases a little bit until about 30,000 years, when it is captured by the 3:2 
mean-motion resonance with the giant planet. It's also the time that the evolution tracks of the 
planetesimal with and without a giant planet in the system diverges. Before the capture point, 
the two tracks almost exactly overlap each other, while after that time, the resonance interaction 
pushes the planetesimal to move outward. Meanwhile, the ratio of the planetesimal's semi-major 
axis to the planet's becomes fixed on some value (in this case, the ratio is 1 .5^). Figure|4]shows the 
evolution of their period ratio. Starting from about 2.3 : 1 the ratio of the two periods decreases 
until the resonance point 3:2, and then it remains fixed on this value. During this process, although 
the period ratio passes 2: 1 at about 10, 000 years, the planetesimal was not captured by the planet's 
2: 1 mean motion resonance. 

As we have mentioned above that there are many factors which determine planetesimal's 
resonant capture probability. First, for planetesimals with certain size, the initially positions of the 
planet and the planetesimals are important. Once the planet's original orbit is set, only a range of 
initial position is appropriate for the planetesimals which would be captured. In Figure [51 we show 
the relation between initial positions and the final positions (semi-major axis) during the orbital 
evolution of several planetesimals with a range of sizes. In this model we adopt the same planetary 
system as in Figure [3l with a linear mass loss prescription for the host star (a constant loss rate 
from Mi = 4 Mq to Mf = 1 Mq within 100,000 yr and Vv„„d = 10 km s"'). A IMj planet is placed 
on an initial circular orbit with a radius 20 AU, and its orbit evolves to rj = 80 AU, according to 
the ratio of the star's mass loss. All planetesimals have density pp = 3 g cm'^. This figure shows 
that planetesimals with size in the range of 3km-30m could be captured onto the planet's mean- 
motion resonance, including 2:1, 3:2 and 4:3 resonant point (period ratio). It shows that despite a 
wide range of initial positions, for planetesimals of the same size, only the ones with initial a, in a 
"critical range" could be captured onto the planet's mean motion resonances. For example, among 
all the s = 300 m planetesimals starting from 25 AU to 37 AU, only the ones with a, between 27 
AU and 31 AU could be captured onto 3:2 mean motion resonance. 
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Second, a planetesimal's size is also a determining factor for resonant capture probability. 
Equation (fT5l) indicates that the effective size range for resonant capture (^rescap < ^ < ^catchup) 
depends on the planetesimal's initial location a,, so the "critical range" of initial positions in the 
first point depends on the size of planetesimals. This dependence is confirmed by our numerical 
calculations. In Figure [5]the "critical range" for 5 = 1 km object to be captured onto 3:2 resonance 
is from 26 AU to 28 AU, while for s = 100 m object is from 26 AU to 38 AU. The smaller 
planetesimals comparatively feel "stronger" gas drag than larger ones at the same orbit, so they 
must have larger initial a,'s so that the gas drag is insufficient to drag them in and they have the 
possibility to being captured by the planet. Generally speaking, this "critical range" is larger and 
further for small particles, as the picture demonstrates. These results are in good agreement with 
the above analytic estimates. 

The magnitude of the gas drag is proportional to Pwmj which is determined by the velocity 
of the gas wind and the mass loss rate of the central star. However, equation © indicate that Sc 
is essentially independent of M* and Vwind- Nevertheless ^rescap does depend on tm and the star to 
planet's mass ratio as well as planetesimals' internal density (through Sc). In order to explore these 
dependence, we consider variation in the model parameters from the standard model we have been 
working with above, i.e. with M ~ 10~'^Mq yr~' and Vvy,,,^ = 10 km s These quantities may vary 
not only from star to star but also temporally around any host star. In addition, the initial and final 
mass of the stars determines the extent of planet's orbital expansion and sweep up of the residual 
planetesimals. The internal density of the planetesimals (pp) may also have a range of possible 
values, from ~ 3gm cm"^ for asteroids to an order of magnitude smaller for comets. 

Apart from the model parameters we have discussed above, we have adopted several different 
prescription for the evolution of the mass loss rate including constant, exponentially increasing, 
or ramp ups followed by declines. We find that once we have specified the initial and final stellar 
mass and the duration of the mass loss, our results are insensitive to variations in the mass loss 
rates. We also consider a range of planetary masses (1 - lOM^) and eccentricities (0 and 0.1). Our 
results are also insensitive to the magnitude of these parameters. 

In the present context, a planet's final semi-major axis is generally less than 100 AU and its 
asymptotic period is < 1000 yr. Thus, their orbital expansion is always adiabatic such that their 
orbits expand by a factor of Mi/Mf. By specifying the planet's initial orbital radius, we can assess 
the inventory of planetesimals to be captured by the planet. 

At which ratio of the orbital circle of the planetesimal and the planet would the particles be 
fixed is also determined by many factors. The rate of the orbital circle of the planetesimal and the 
planet would keep rising while no resonance happens. When it close to a resonance point, whether 
or not it could stop raising and stay is primarily determined by whether or not the resonance is 
strong enough, which is decided by factors like the gas' density, the positions of the planetesimal 
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and the planet, etc. In our discussion, there are three major resonance points: rate of the orbital 
circle equals to 3:2, 4:3 and 3:2, while 3:2 is the most important one. It's shown in Figure[5l Take 
the curve of s = 300m planetesimals for example: the objects starting from a, =27 AU-31 AU are 
captured onto 3:2 resonance, with ay around 105 AU; starting from a,- =25 AU-26 AU are captured 
onto 4:3 resonance, which have a/ around 97 AU; and starting from a, =35 AU-36 AU are captured 
onto 2:1 resonance, which have a/ around 128 AU. 

5. CoUisional evolution of planetesimals and production of grains 

In the previous section, we suggest that km-size icy particles are the preferred population 
of residual planetesimals for planets' resonance capture during the AGB mass loss phase of their 
host stars. In the next section, we first estimate the condition under which a population of such 
planetesimals may emerge at 10-100 AU prior to the AGB phase. 

We then show that as a consequence of their orbital expansion and planets' sweeping pro- 
cess, these planetesimals' eccentricities increases to modest values and their collision frequency 
is strongly enhanced. We show that planetesimals orbit crossing with each other leads to frequent 
high-velocity disruptive impacts and the production of small grains. These grains can then repro- 
cess the stellar radiation and produce the observed FIR rings around the white dwarf at the center 
of the helix nebula. 



5.1. Inventory of planetesimals 

In §2, we indicated that around typical progenitor of planetary nebula (A main sequence star), 
there may be a rich population of planetesimals. We now estimate their characteristic size prior to 
the mass loss from their host star. 

There are two possible formation epochs for residual planetesimals (first few Myr and through 
out the main sequence life time of their host stars r*). First-generation planetesimals are formed in 
gas-rich protoplanetary disks. In such a gaseous environment, the velocity dispersion of the small 
planetesimals (cr) is excited by their mutual interaction and damped by their interaction with the 
disk gas (Kokubo & Ida 2002). Due to hydrodynamic drag by the turbulent disk gas, these particles 
also undergo orbital decay (Supulver & Lin 2000). During this formative stage for their host star, 
largest planetesimals can grow to sizes well beyond a few km, especially in special locations such 
as the snow line and in the proximity of proto gas giants. They can also be scattered to large 
distances by the emerging planets (Zhou & Lin 2007). Nevertheless, we have already showed in 
the last section that during the planetary nebula stage, planetesimals with s > ^'catchup retain their 
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semi-major axis ratio with respect to any embedded planet. They are unlikely to be captured onto 
any planet's expanding mean-motion resonances. 

Planetesimals can also grow in situ during the main sequence life span of their host stars 
(Th. ~ 10^"^ yr), long after the depletion of the disk gas. In this case, their cr is determined by their 
own coUisional properties. Collisions can lead to both mergers and fragmentation. The condition 
for coUisional fragmentation for basalt and ices have been extensively studied experimentally and 
simulated numerically (Stewart & Leinhardt 2008). In collisions between two comparable km-size 
aggregates, the critical impact speed Vfrag for significant fragmentation is around a few meter per 
second (which is comparable to their surface escape speed Vg for these planetesimals). Shattering 
collisions are not only inhibits the coagulation of planetesimals but also introduces considerable 
energy dissipation. In Keplerian disks (such as Saturn's rings), particles' velocity dispersion (cr) 
is determined by its excitation due to elastic scattering and damping due to inelastic collisions 
(Goldreich & Tremaine 1977, Bridges et al 1984). In gas-free disks with cr < Vfrag, planetesimals 
undergo cohesive collisions with growing a (Palmer et al. 1993, Aarseth et al. 1993). When their 
cr exceeds Vfrag, collisions quench both mass and cr growth. In the absence of a dominant popula- 
tion of large planetesimals or embedded planets, planetesimals establish a coUisional equilibrium 
in which their cr is approximately isotropic with cr ~ v^ag- An isotropic velocity dispersion implies 
that planetesimals can only collide with others in an annulus with a half width ajCr/vK and the mag- 
nitude of cr determines not only the coUisional frequency but also the thickness of the planetesimal 
disk. 

For km-size planetesimals, v^ag ~ such that the coUisional cross section of these planetesi- 
mals is mostly determined by their physical size. If most collisions are cohesive, the growth time 
scale (Tgrowth) would be comparable to the coUisional time scale (Tcoi) 

Tgrowth - Teoi ^ (aj sfiMJ M,,d(P K /2n) (18) 

where and Mtot are the individual and characteristic total (at a,) masses of planetesimals. The 
present estimates for Mtot in the Kuiper Belt ranges from a fraction to a few and that of debris 
disks around A stars may be near the upper end of this range (see §2). From equation (fT8l) , we find 
the characteristic size of planetesimals to be 



The above equation suggests that the spontaneous production of an isolated population of km-size 
planetesimals within the main sequence life span of a A star requires Mtot to be a few M^. This 
estimate is consistent with that simulated for the Kuiper Belt region (Kenyon & Bromley 2004). 

In the derivation of 5char, we assume that Tgrowth ~ Tcoi, i.e. all collisions are cohesive. Since 
o" ~ Vfrag some collisions are disruptive and ^chai represents an upper limit in the planetesimals' 
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size distribution. In §2, we briefly show that the particle size distribution may follow a power law 
in which 

dNIds = (A^o/Vax)(^/Vax)"'' (20) 

where ^^ax is the size of the largest particles and A^^ is a normalization factor. With this size 
distribution, most of the mass are attained by the largest planetesimals. As long as 5char > 1 km> 
there is always a supply of residual planetesimals which may be captured by the resonance of 
expanding planets. 

Equation (fT9l ) implies a large potential dispersion in 5char around stars with diff"erent r* and a,. 
Perhaps the largest uncertainty is the total mass of residual heavy element Mjot. Around stars with 
relatively massive debris disks, Schm > 1 km and the total mass of km-size planetesimals is 

Mkn, - (lkm/w)'^'Mto, M/„f . (21) 

But, in less massive disks, the supply of km-size planetesimals may be limited. Based on this 
consideration and the uncertain fraction of A stars with planets, we do not expect the detection of 
trapped planetesimals around the remnant white dwarf at the center of every planetary nebula. 



5.2. Eccentricity excitation 

After km-size planetesimals are captured onto the planet's mean-motion resonances, their 
orbital expansion would be locked provided the orbital migration time scale is longer than the 
resonant libration time scale i.e. tm > Tub or equivalently, 

Mp > (M,/4expl)(/'^/rM)' - IQ-'M,. (22) 

Similar processes are found in the tidal evolution of Galilean satellites (Goldreich 1965, Lin & 
Papaloizou 1979, Peale & Yoder 1981) and migration of resonant planets (Lee & Peale 2002, Kley 
et al. 2004). 

During their forced co-moving orbital expansion, the conservation of an adiabatic invariance 
induces the growth of resonant planetesimals' eccentricity (C s) (Murray & Dermott 1999) on a time 
scale 

Te,e = eje, = 2ejTexp = 2e]TM- (23) 

The above expression is for planetesimals in the planet's 2:1 mean motion resonance. A similar 
expression can be derived for planet's other mean-motion resonances. 

Although the stellar outflow is in the radial direction, its drag on the planetesimals can also 
lead to the damping of their eccentricities on a time scale 

Te,d - is/Sc)TM. (24) 
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The above expression is derived from a linear perturbation analysis of equation ^ in the limit of 
small e^- Comparison between Te,e and ^, it is apparent that planetesimals' grows initially until 
an asymptotic equilibrium is established in which 

ee - {sls,fl\ (25) 

After they are captured by the planet at a location Vp = reap, the planetesimals' eccentricity growth 
may also be quenched by the termination of the stellar mass loss with an asymptotic 

ea - (ln(rfi„,i/reap))^/'. (26) 

In general, the eccentricity of co-moving resonant planetesimals is excited to the minimum 
of Ce and which may both be a significant fraction of unity. (Both of these estimates on the 
asymptotic eccentricities break down as their values approach to unity.) Since the semi-major axes 
of the resonant planets are in lock-steps with that of the planet, the orbits of widely separated 
planetesimals overlap. The onset of orbit crossing creates opportunities for collisions between 
these planetesimals. 

The eccentricities of distant planetesimals with period Pk ~ tm (at a, > 500 AU) may also 
increase as a consequence of their marginally non-adiabatic adjustment to the host star's declining 
gravitational potential. However, these planetesimals' peri-center distances would not decrease 
significantly so that they are unlikely to venture into the regions of interest. At these large distances, 
planetesimals are sparsely populated, their collision frequency is also expected to be low. 

In order to verify these analytic estimates, we show in Figure [6] the evolution of the planetes- 
imal's eccentricity in the system in Figure [31 In this case, it remains negligible until 10,000 yr. 
when it encounters the planet's 2:1 mean-motion resonant point. Although not being captured, 
its eccentricity still slightly increases. When the planetesimal is captured onto the planet's 3:2 
mean-motion resonance at 30, 000 years, its eccentricity is evidently excited. At the end of our 
simulation it reaches its asymptotic value ~ 0.24. These result is in good agreement with the above 
discussions. 



5.3. CoUisional frequency and outcome 

The results in Figure [6] can be generalized. In Figure |7l we show the asymptotic eccentricity 
distribution of a population of planetesimals (the system in Figure [3]) verse their final orbital semi- 
major axis. It's very clear, as we pointed out before, there are three major resonance points of this 
system which can capture planetesimals onto co-moving orbits with giant planet and excite their 
eccentricities: 3:2 resonance point (the most important one, final semi-major axis of planetesimals 
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ttf ~ 105 AU); the 2:1 resonance point (a/ ~ 97 AU) and 4:3 resonance point (a/ ~ 128AU). These 
three resonance points correspond to the three vertical bands on the plot. Once being excited, the 
eccentricity of the planetesimals will increase to about 0.1 ~ 0.3. For planetesimals which are 
captured onto 3:2 resonant orbit, since their asymptotic semi-major axis is around 100 AU their 
orbits cross each other in the region between ~ 80 - 120 AU. Thus the planetesimals' eccentricity 
excitation enlarges the radial extent of their impact zone from ~ a,Vfrag/v/f ~ 10"^a/ to ~ Icsai ~ 

0. 4a,. This expansion greatly increases the pool of potential impactors. 

As planetesimals along the expanding paths of the planet are captured onto its mean-motion 
resonances and co-moved by it to large distances, the surface number density 2p of planetesimals 
can increase by a factor (~ a few), depending on the initial distribution of prior to the star's 
mass loss. (In this estimate, we took into account the growth of planetesimals' eccentricity associ- 
ated with their resonantly forced orbital expansion.) 

In contrast, the planetesimals' dispersion velocity in the direction normal to the disk plane 
does not change significantly (from its initial value cr ~ Vfrag) because the forcing by the mean- 
motion resonances is applied in the radial direction. (The velocity dispersion is no longer isotropic, 

1. e. the radial component of the velocity dispersion ellipsoid, cr^ ~ e^VK, becomes much larger than 
that in the direction normal to the orbital plane, <t, ~ Vfrag-) Consequently, the spacial density of 
resonant planetesimals {n'^) actually increases from their initial values {Mtot/(Me2najH) where H ~ 
cr^ai/vK is the disk thickness) by a factor similar to fz and the coUisional frequency is enhanced 
by both the density and velocity enhancement. The characteristic coUisional time scale Tcoi for the 
km-size planetesimals (which are trapped in the planet's mean-motion resonances) is reduced from 
its original magnitude to 

<oi - U(n'pO-rS^) - TUvfr.Je,VKk){l s,)"^ - (10-^ - 10-^)(1 km/5,)i/2^,. (27) 

In the above expression, we assume that within the main sequence life span (t*) of the host star,z.e. 
prior to the AGB mass loss phase, perfectly cohesive collisions led to the formation of planetesi- 
mals with the power-law size distribution in equation (l20l) up to a size Sc- 

During and after the AGB mass loss phase, the radial velocity dispersion associated with the 
asymptotic mean eccentricity (~ 0.2) is ~ 0.5km s~' at ~ 100 AU. Since this value of a is much 
larger than v^ag, we anticipate all collisions during this stage will lead to powderization of the 
planetesimals. Since only km-size planetesimals are resonantly captured and co-move with the 
orbital expansion of the planet, production rate of fragment would be 

Mfrag - (1 km/5,)'/2^tot/T;„i - (e.vj,A/vfrag)(Mtot/T.) ~ lO^-^Mtot/r,). (28) 

With this erasion rate, the total mass of fragments produced during the AGB mass loss phase 
would be 

M[,, - MfragTM - MU^Q^-'tmIt,). (29) 
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During the mass loss epoch, the magnitude of r^^j can be reduced to ~ 10"*"^ yr which would be 
comparable to tm- In these systems, a significant fraction of Mtot could be fragmented through 
these disruptive processes. Since the emergence of km-size planetesimals requires Mtot to be a few 
Me, we anticipate that in planetary nebula where they coexist with some gas giants, the total mass 
of their fragmentary debris (Mt'^t) may indeed exceed that of the Earth. 

Fragments' mass spectrum, immediately after the disruptive impacts, depends on the highly 
uncertain material strength of planetesimals. If these planetesimals are loosely bound aggregates 
high speed (~ km s"^) impacts may disperse them down to the sizes of their basic building blocks. 
For simplicity, we assume a power-law size distribution as expressed in equation (l20l) . In this 
case, a major portion of M[^^ is contained in the largest fragments (with s = Smax and mass 
'"max = 4;7rpp5^3x/3) SO that the normalization factor can be derived from equation (|20l ) to be 
No = M[^J2m^^^. In the absence of reliable experimental data, we infer below the magnitude of 
5max using the flux of reprocessed radiation from the ring around WD 2226-210. 



5.4. Reprocessed radiation by the residual fragments 



Another implication of the power-law size distribution is the fraction of stellar photons which 
may be absorbed and scattered by the fragments. In the calculation on the intensity of reprocessed 
radiation at any instance of time, we need to take into consideration the potential depletion of all 
particles. We have already indicated above that, during the gas loss phase, fragments with sizes 
s < ■^blowout are continually entrained and removed by the outflowing gas from the location of their 
production on a time scale Tfly ~ ai/Vwind- Similarly, inward orbital decay of fragments in the 
size range i'biowout < s < ^bypass proceeds on a time scale Tdrag = (s/Sc)tm (equation (flOl) ) are not 
inhibited by the planet's tidal barrier during their orbital decay. 

In the calculation of intensity of reprocessed stellar irradiation at a particular wavelength range 
(A), we also need to include a reduction factor (s/A) for emissivity of small particles (in the range 
from = A> s down to a minimum grain size i'min). In the limit ^bypass ~ > ■^blowout > ■^.i > ■^min, 
the total area covered by the fragments in the impact zone would be 
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Neglecting the YORP eff'ect and inserting expressions for dN/ds in equation (l20l) . we find that 
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Although AGB stars' L* can reach ~ 10 '^Lq to exert non-negligible YORP effect on km-size 
planetesimals at ~ 10-100 AU, this phase is brief and does not significantly affect their orbits. 

In the limit Sc ~ ^bypass, CTzone ~ ^^^''■5'ma^x ■^bypass- Interior to the impact zone and the planet's 

orbit, fragments with ^bypass > ^blowout > Sa contributes to a-interior ~ (27rA^>^ax/*c)[(*bypass/*max)'^^- 

(^biowout/^max)'''^] ~ 2nN'„slil,sll^'^^^^ < (T^one bccausc the normalization factor A^^ (which is deter- 
mined by the resonant capture efficiency) is expected to be smaller than No- Exterior to the impact 
zone, small fragments (s < ^blowout) contributes to cTexterior ~ (27rA^o5Lx'^fly/TM)[2(Vax/*i)^^^ - 
(■^max/ ■^blowout) (.^max/ ^a)(. Smm Vax)^''^] ~ o-zone(*bypass/*^)(Tfly /tm) wMch may also bc much 
smaller than (Tzone- 

Equations (fTTl) and (fT6l) indicate that both i'biowout and ^bypass increases with M,. After the host 
star attains its Mf (as in the case of the white dwarf WD 2226-210 at the center of the Helix neb- 
ula), the mass flux of gas outflow is significantly reduced. While the eccentric resonant-captured 
planetesimals continue to collide with each other on a time scale r^^j and generate fragments at a 
rate Mfrag- When both i'biowout and ^bypass decreases below sa, equation (|3T1) is much simplified by a 
continuous power-law size distribution (as in equation [201) to 



C"zone ~ ^TtA^o i'jjjax 



1/2 



1 



In this case, the fragments are confined to the impact zone. 



(32) 



Finally, the particle size distribution may also be truncated at the lower ranges the radiative 
blow out process (Sy). In that limit, would be replaced by Sy in the above expressions. 



6. Structure of the disk around WD2226-210 

The proximity of WD2226-210 at the center of the Helix nebula has been imaged in multiple 
wavelength. The observed spectral energy distribution (SED) has been analyzed in terms of a disk 
of dust grains which are assumed to have the same astronomical size distribution as equation ( |20l ) 
and a size range between 60-1000 yum (Su et al. 2007). (The lower size limit is set (^min = 60/um) 
by the condition for radiative blow out (i.e. Sy), albeit ^^ax is arbitrarily set.) The total area covered 
by these dust grains is inferred to be crtot ~ 6x10^^ cm^ distributed over a distance between 35-150 
AU. For the assumed range of grain sizes, their total mass is inferred to be Mgrain ~ O.lMe. (This 
mass is inferred from the assumed size distribution in equation [201 neglecting the contribution of 
hydrodynamic blow out and drag on the small grains in equation [3T|). The magnitude of 

Mgrain - 0.1(Vax/0.1cm)i/2Me (33) 
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if the assumed size distribution extends to grains larger than 1000 ;um=0.1 cm.) Since the present 

and past PR drag time scale for these grains (~ 10 My) is < r,, they must be continually replen- 
ished, presumably through the powderization of their parent-body residual planetesimals. Su et al. 
(2007) estimated the current collision time scale for lOO/im grains to be ~ 0.1 - 1 My which is 
> tm. 

Although this model provides a reasonable fit to the observed SED, a challenging theoretical 
issue is whether the grains' parent bodies may have sufficient collision frequencies to generate so 
many fragments despite their (parent bodies') much smaller surface area to volume ratio. If the 
generation of Mgrain requires a substantially more massive (>> M®) population of parent bodies, 
its protracted retention during the main sequence life span of their host stars would need to be 
addressed rather than assumed. 

We now apply the results in the previous sections to provide some quantitative justification 
for the model proposed by Su et al. (2007). In §2, we cite evidences that both gas giant planets 
and debris disks are commonly found around A stars. In our model, we consider the possibility 
that the progenitor star of WD 2226-210 was surrounded by one or more gas giant planets (at 
distances comparable to those around Uranus and Neptune) and a debris disk similar to the Kuiper 
Belt in the solar system. In previous chapters, we have already showed that in debris disks which 
contain a few (or more) M® within ~ 100 AU, a population of km-size planetesimals (with a total 
mass Mrmtot > 1 Md,) can form through inelastic and cohesive collisions (with a velocity dispersion 
c ~ Vfrag ~ a few m s~') during the main sequence life span of type A stars (i.e. within ~ 100 
My). Due to the combined planetary tidal perturbation and the YORP eff"ect, a fraction of these 
residual planetesimals may accumulate near the planets' outer mean motion resonances. 

During the AGB mass loss phase when most of the stellar mass is rapidly lost, semi-major 
axes of both planets and relatively large residual planetesimals increases as M^'. In contrast, the 
outflowing wind imposes hydrodynamic drag on the planetesimals and causes those with relatively 
small sizes (s < Sc ~ 1km) to undergo orbital decay. When they cross the planet's expanding path, 
these small planetesimals are captured onto its mean motion resonances of the planets. 

During the subsequent shepherd orbital expansion, the resonant planetesimals are swept up 
with enhance surface densities (by a factor /s ~ a few), excited eccentricities (e^ ~ or 
which can be up to a significant fraction of unity), and highly anisotropic velocity dispersion 
(cTr ~ CsVk ~ 10^"^ (T^). The coUisional time scale for these planetesimals are greatly reduced (by 
a factor ~ 10^""^) from t* to ~ tm so that a significant fraction of the planetesimals captured on the 
planet's mean-motion resonances may actually collide with each other. The impact speed of these 
collisions (~ cr^. ~ 0.5km s~^) is sufficiently large to pulverize km-size planetesimals to very small 
pieces. 
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The actual size of the largest fragments (i'max) is poorly determined. However, we can infer it 
by assuming Mgy^m - ^tot (from equations [33] and [291) and (Xtot ~ cr^one (from equations [311 or [321) . 
If we neglect any residual stellar winds and self-consistently adopt M^^^ ~ Mjot = 1 M®, we would 
infer Sjnax ~ 10cm. 

In the absence of significant stellar wind, hydrodynamic blowout and drag are negligible. 
Whereas radiation blow out remove smallest grains (with s < 60jum), the PR drag causes intermediate- 
size (~ 100//m) grains to undergo orbit decay on a time scale (tpr) of 3-50 Myr over this region 
(between 35-150AU). Since this time scale is much longer than the time lapse since the AGB mass 
loss phase around WD 2226-210, it has negligible effect on the ring at the center of the Helix 
nebula. But around mature white dwarfs (with r* > tpr), it is possible for 

Ta(A = Rr, Tdecay = ^Pr) > Tsyn (34) 

or equivalently, tpr > {M^IMpfl^Px- In this limit, these grains may not be able to bypass the 
planet's orbit so that their asymptotic surface density distribution (around older white dwarfs) may 
bear the signature of planet's tidal barrier and have an inner edge (Wyatt et al. 1999). Nevertheless, 
one or more eccentric planets can disable the orbits of these planetesimals, especially after a sub- 
stantial loss of (Levison et al. 1994) This process may also provide a supply of heavy elements 
to pollute the atmosphere of DZ white dwarfs (Debes & Sigurdsson 2002). We will consider that 
effect elsewhere. 

In the case of WD 2226-210 at the center of the Helix nebula, there may still be a reduced 
but non negligible stellar wind. Due to hydrodynamic drag (rather than stellar photons), sub-mm 
grains are may be blow away while slightly larger grains may undergo orbital decay and bypass any 
embedded planets. In this limit, grain depletion would lead to transitions in their size distribution. 
The inferred values of imax could be reduced by the depletion of the small particles. 

Finally, we consider possible observable consequence of this process. During the planetary 
nebula, embedded planets may also accreted gas from the outflowing stellar wind at a rate 

M,.(^f(^VM.. ,35) 

The above equation is derived from a modified Bondi accretion formula. For a Jupiter mass planet, 
the total accreted gas during the AGB mass loss phase (~ M^tm) would be a small fraction (~ 
Mp/Mt) of its own mass. Nevertheless, it could provide an accretion luminosity (up to ~ 10~'^Lq) 
which may be marginally detectable by ALMA against the ring background. 

The hot white dwarf WD 2226-210 at the center of the Helix nebula provides an intense 
source of UV radiation. The ionization of planetary atmosphere can lead to its photoevaporation 
and mass loss (Lammer et al. 2003, Baraffe 2004, 2005, Murray-Clay et al. 2009). Lyman 
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a absorbing clouds have been observed around a close-in planet, HD 209458b (Vidal-Madjar et 
al. 2003). Its extension beyond the planet's roche lobe has been attributed to be a signature 
of outflow, albeit the inferred mass loss rate remains uncertain (Ben-Jaff'el 2007). The energy 
flux of ionizing photons irradiated onto the surface of gas giants located at ~ 30 AU from WD 
2226-210 is 2-3 orders of magnitude larger than that onto hot Jupiters around solar type stars. 
Scaling from the radiative hydrodynamic models for HD 209458b (Murray-Clay et al. 2009), we 
do not expect long-period gas giants to loss significant fraction of their mass during the AGB mass 
loss phase. Nevertheless, photoionizing stellar photons may induce emission of signature spectral 
lines on the extended planetary envelope which may potentially be detectable with high dispersion 
spectroscopic observations (we thank Y.-H. Chu for this suggestion). 

Mean motion resonant capture of km-size planetesimals by one or more embedded planets 
and their subsequent lock-step orbital migration also lead to their non axisymmetric surface density 
distribution (Murray-Clay & Chiang 2005). In principle, such a structure can be resolved with sub 
arc second imaging. However, the observed IR radiation is due to reprocessed stellar radiation by 
the fragments of high-speed collisions. Their recoil velocity may be sufficient to smear out such 
a structure. However, multi- wavelength resolved images may reveal a particle size dependence in 
the grains radial distribution. The detection of very small is < .Vbiowout) grains at large disk radii 
and intermediate s (^blowout < * < ^bypass) would provide supporting evidence for our scenario. 



7. Summary and discussion 

In this paper we construct a dynamic model for the post main sequence evolution of planetary 
systems. With this model, we find that as a consequence of mass loss from the central star, km- 
size planetesimals can be captured by mean-motion resonance and form a disk beyond the planets. 
During the subsequent orbital expansion, the eccentricity of the resonant planetesimals is excited 
which enhances their collision frequency. Highly destructive impacts among these planetesimals 
can continuously provide a ring of dust around the mature star, which can be seen in the FIR and 
MIR wave band. This phenomenon can also be extended to other A-F stars, with the inference 
that a large fraction of stars may have planet systems and substantial debris disks during their main 
sequence epoch. Stars with higher main sequence mass, greater mass loss rate and lower velocity 
of star wind tend to have higher ability to capture the planetesimals. 

We applied this model to account for the dust ring around WD 2226-210 at the center of the 
Helix nebula. It can also be used to study other young white dwarfs. For example, we anticipate this 
process will occur when the Sun evolve to its AGB mass loss phase, albeit among of fragmentary 
dust would be limited and the intensity of the reprocessed radiation may be weak. Since this 
process depends sensitively on the initial heavy element content and the presence of sufficiently 
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massive planets, we expect there to be large dispersion in the signatures of dust debris around 
young white dwarfs. Nevertheless, we expect to see more dusty rings around planetary nebulae 
with A star progenitors because the fractions of stars with planets and debris disks is observed to 
increase with stellar mass. 

The present work is develpoed for the pwind °^ 1 /r^ (or v^ind = constant) portion of a stellar 
wind and thus not applicable to the region very close to the star. In addition, apart from the plan- 
etesimals that are captured by the planet, the fate of other planetesimals also need to be discussed. 
These may also play an important role in the evolution of the white dwarf. 

Some observational tests of these models should be possible with present technology or will 
be possible soon. As the planetesimals are collected by the giant planet, it is most likely that the 
gas giant planet survives the AGB period, and the orbit changes to several times the original one. 
If such giant planet around white dwarf can be detected, it would provide support to our model. 
Although there are several failed attempts in these searches (Debes et al. 2005a, b, 2007), the 
upper mass limits set for the unseen planet is far larger than that needed in our model. 

We thank Sverre Aarseth for provide the Hermit integrator, Y.-H. Chu and an anonymous ref- 
eree for useful correspondence and helpful suggestions. This work is supported by NSFC(10233020), 
NCET (04-0468), NASA (NNX07A-L13G, NNX07AI88G, NNX08AM84G), JPL (1270927), and 
NSF(AST-0908807). 
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Fig. 1. — Orbital evolution of two planetesimals around a mass-losing star. The initial mass of 
the central star in this simulation is 4 Mq. It loses 75% of its mass at a constant mass loss rate 
over 10^ years (from now on it is called "linear mass loss"). The upper and lower curves in the 
plot show the evolution of semi-major axis verse time of two identical planetesimal, with density 
p = 3 g/cm^ and size 5=1 km. The upper one starts from 15 AU, and its size is larger than the 
Sc at that position, so its orbit expands. While the lower one starts from 10 AU, with its size being 
smaller than Sc at that position, it undergoes orbital decay, and finally drops into the star within 
about 25,000 years. 
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Fig. 2. — Orbital evolution of 4 planetesimals. The mass of the star changes linearly from an 
initial value M, = 3 Mq to a final value Mf = 0.5 Mq in 50,000 years. Curves show the evolution 
of semi-major axis as a function of time for 4 objects all starting from 15 AU, with size s = (from 
top to bottom) 10 km, 1 km, 100 m and 10 m. The wind velocity is v„ind =10 km s"' 
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Fig. 3. — This plot compares the orbital evolution of a * =0.1 km planetesimal when there is a 
giant planet in the system (red) or not (green), while the blue curve shows the orbital evolution of 
the planet. The density of the planetesimal is = 3 g cm"^ and it starts out with a circular orbit 
at 35 AU around a M, = 4 A/q star which linearly reduces its mass to Mf = 1 Mq on a time scale 
of 100,000 yr through a wind with a speed Vmnd = 10 km s'K The planet with a mass Mp = I Mj 
was placed in a circular orbit at 20 AU initially. 
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Fig. 5. — The orbital evolution of a group of different-size planetesimals verse time. The size of 
the planetesimals in the five curves are (from left to right) 3 km, 1 km, 300 m, 100 m and 30 m. 
The system is the one in Figure [31 
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Fig. 6. — The eccentricity evolution of the lOOm-size planetesimal in the system of Figure[3] 
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Fig. 7. — The final eccentricity verse final semi-major axis for several planetesimals with diff'erent 
sizes in the system of Figure [51 



